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Abstract. We provide a rigorous treatment of the inverse scattering trans- 
form for the entire Toda hierarchy for solutions which are asymptotically close 
to (in general) different finite-gap solutions as n — > ±00. 



1. Introduction 

The Toda lattice is one of the most prominent discrete integrable wave equations. 
In particular, it can be solved via the inverse scattering method. For the classical 
case, where the solution is asymptotically equal to the (same) constant solution, 
this is of course well understood and covered in several monographs (e.g.) [T^ . 
[36j . or [33] . The corresponding long-time asymptotics were first computed by 
Novokshenov and HabibuUin [29 and were later made rigorous by Kamvissis [16] 
under the additional assumption that no solitons are present (the case of solitons 
was recently added in [23]; see also the review f22]). 

The inverse scattering transform for the entire Toda hierarchy in the case of a 
finite-gap background was solved only recently by us in [5] as a continuation of 
[31j . Similar results were obtained by Khanmamedov [21j . Long-time asymptotics 
for such solutions have been given by Kamvissis and Teschl for the case without 
solitons [H], [H], [20] and by Kriiger and Teschl for the case with solitons [24| (for 
related trace formulas and conserved quantities see [27]). 

In this respect it is important to mention that even the important case of a 
one-soliton solution on a finite-gap background has different spatial asymptotics as 
n ±00 and hence is not covered by the above results (see [11], [34] )• Hence this 
clearly raises the need to extend the results from to the case of solutions which 
are asymptotically equal to (in general) different finite-gap solutions as n ^ ±00. 

In fact, the simplest case, where the solution is asymptotically equal to two 
different constant solutions, has already attracted considerable interest in the past. 
The first to solve the corresponding Cauchy problem (in the case of rapid decay 
with respect to the background) seems to be Oba [30]. Moreover, the long-time 
asymptotics were considered in [1] , [3] , [4] , [14] , [15] , [T^ , . 

Our aim here is to fill this gap and to provide a treatment of the inverse scattering 
transform for the entire Toda hierarchy in the case of steplike quasi-periodic finite- 
gap backgrounds. Note that since we treat the entire Toda hierarchy, our results 
also cover the Kac-van Moerbeke hierarchy as a special case [28] . 
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Finally, wc remark that the corresponding result for the Korteweg-de Vries equa- 
tion is much more involved and was only recently solved by Egorova, Grunert, and 
Teschl ,6 under some additional restrictions on the spectra of the background op- 
erators. 

After introducing the Toda hierarchy in Section [2l we will first show that a 
solution will stay close to a given background solution in Section [3l This result 
implies that a short-range perturbation of a steplike finite-gap solution will stay 
short-range for all time, and it shows that the time-dependent scattering data satisfy 
the hypothesis necessary for the Gel'fand-Levitan-Marchenko theory ^25^. This 
result constitutes the main technical ingredient for the inverse scattering transform. 
In Section |4] we review some necessary facts on quasi-periodic finite-gap solutions 
and in Section[5]we compute the time dependence of the scattering data and discuss 
its dynamics. 

2. The Toda hierarchy 

In this section we introduce the Toda hierarchy using the standard Lax formalism 
([25]). We first review some basic facts from [2 (see also [13], |33j). 
We will only consider bounded solutions and hence require 

Hypothesis H.2.1. Suppose a{t), h{t) satisfy 

ait) e £°°(Z, M), b{t) e £°°{Z, M), a{n, t) ^ 0, (n, t) G Z x M, 

and let t h-^ {a{t),b{t)) be differentiable m 1°°{'L) ^°°(Z). 

Associated with a{t), b{t) is a Jacobi operator 

(2.1) H{t):e {1)^1^1.), f^T{t)f, 
where 

(2.2) Tit) fin) = ain, i)/(n + 1) + a(n - 1, t)/(n - 1) + bin, i)/(n) 

and (Z) denotes the Hilbert space of square summable (complex- valued) sequences 
over Z. Moreover, choose constants cq = 1, Cj, 1 < j < r, c^+i = 0, set 

j 

gjin,t) = y^^Cj-e{Sn,Hit)''dn), 

3 

(2.3) /ij(n,i) = 2a(n,i)^c,_f(5„+i,ff(<)^5„) +c,+i, 

f=0 

where ((5m,^<5n) denote the matrix elements of an operator A with respect to the 
standard basis, and consider the Lax operator 

r 

(2.4) P2,+2(i) = -HitY^^ ^(2a(%(t)5+ - h,it))HitY-^ ^ g^+xit), 

3=0 

where 5^/(n) — /(n± 1). Restricting to the two-dimensional nuUspace 

Ker(r(t) - z), ze C, 
of rit) — z, we have the following representation of P2r+2(0- 



(2.5) 



P2r+2it) 



= 2ait)Griz,.,t)S+ -Hr+iiz,.,t), 

Ker(r(t)-z) 
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where Gr{z, n, t) and Hr+i{z, n, t) are nionic polynomials in z of the type 

r 

Gr{z,n,t) ^^z^gr-j{n,t), 

r 

(2.6) Hr+i{z, n, t) — z^^^ + z^hr-j{n, t) — gr+i{n, t). 

A straightforward computation shows that the Lax equation 

(2.7) j^H{t) - [P2r+2{t), H{t)] =0, t e M, 
is equivalent to 

TLr{a{t),b(t))i = a{n,t) ~ a{n,t)(^gr+i{n + l,t) - gr+i{n,t)j = 0, 

(2.8) TLr{a{t),b{t))2 = b{n,t)~ (hr+i{n,t)-hr+i{n-l,t)^ =0, 

where the dot denotes a derivative with respect to t. Varying r S No yields the 
Toda hierarchy TLr(a,5) — (TLr(a, 6)i, TLr(a, 6)2) — 0. We will always consider r 
as a fixed, but arbitrary, value. 

Finally, we recall that the Lax equation (j2.7p implies existence of a unitary 
propagator Ur{t,s) such that the family of operators H{t), t e M, are unitarily 
equivalent, H{t) = Urit, s)H{s)Uris,t). 

3. The initial value problem 

First of all we recall the basic existence and uniqueness theorem for the Toda 
hierarchy (see, e.g., [31], [32], or [331 Section 12.2]). 

Theorem 3.1. Suppose (ao,6o) G M = £°^{Z) £°°(Z). Then there exists a 
unique integral curve t 1— > (a(t),b{t)) in C°°(R, M) of the Toda hierarchy, that is, 
TLria{t),b{t)) = 0, such that (a(0),6(0)) = (ao,6o). 

In [31] it was shown that solutions which are asymptotically close to the constant 
solution at the initial time stay close for all time. Our first aim is to extend this 
result to include perturbations of quasi-periodic finite-gap solutions. In fact, we 
will even be a bit more general. Set 

i/p 



(3.1) \\ia,b)\Up={\n 



, l<p<oo, 

(n)(|a(n)| + |6(r.)|). 



sup w[n) [ \a[n)\ + \b(n)\ ] , p — 00 

Then 

Lemma 3.2. Let w{n) > 1 be some weight with sup„ ( | -^^^^j^^p- 1 -I- 1 J^^^i^ I) < 00 and 
fix some 1 < p < 00. Suppose a{n,t), b{n,t) and a±{n,t), b±(n,t) are arbitrary 
bounded solutions of the Toda hierarchy and abbreviate 

(3.2) a(n,t) = " ^ °' 6(n, t) = " " °' 

\a^{n,t), n < 0, I6_(n,t), n < 0. 

Then, if 

(3.3) II (a(i) - d{t),b{t) ~ b{t)) |U,p < 00 
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holds for one t = <o G M, then it holds for all t CzR. 

Proof. Without loss of generality we assume that to = 0. Let us consider the 
differential equation for the differences S{n,t) — (^a{n,t) — a{n,t),b{n,t) — b{n,t)^ 
in the Banach space of pairs of bounded sequences S = (^2) for which the norm 
||<5||u),p is finite. We claim that 6 satisfies an inhomogeneous linear differential 
equation of the form 

S{t)= Ar,jit)iS+yS{t)+Br{t) 

\j\<r+l 

(see e.g. [5] for the theory of ordinary differential equations in Banach spaces). Here 
S^{6i{n, t),S2{n, t)) = {Si{n±l,t),S2{n±l,t)) are the shift operators, 

111 j-\ 112 



^rA^^t) - yAl^^^n,t) Ali{n,t)) 



are multiplication operators with bounded two by two matrix-valued sequences, 
and 

is a vector in our Banach space with Br^i{n,t) = for \n\ > [|J + 1. All entries of 
Arj{t) and Br{t) are polynomials with respect to {a{n + j, t), b{n + j, t)), {a±{n + 
j,t),b±{n + j,t)), \j\ < [|J + 1. Moreover, by our assumption the shift operators 
are continuous, 

P e [1,00), 




and same is true for the multiplication operators Arj (t) whose norms depend only 
on the supremum of the entries by Holder's inequality, that is, on the sup norms 
of {a{t),b{t)) and {a±{t),b±{t)). Finally, recall that by unitary equivalence of the 
operator family H{t), respectively H±{t), we have a uniform bound of the sup 
normsup.^{\a{n,t)\ + \b{n,t)\) < 2\\H{t)\\ = 2||ff(0)||, respectively sup„(|a±(n, t)| + 
\b±{n,t)\) < 2\\H±{t)\\ — 2||iJ-|-(0)||. Consequently, there is a constant such that 
S|j|<r+i ll^'-j'lOII IK'^^)"' II — ^r- Moreover, we will show below that the vector 
Br{t) has only finitely many nonzero entries and thus ||i?r(Ollii',p ^ where the 
constant again depends only on the sup norms of {a{t),b{t)) and {a±{t),b±{t)). 
Hence 



\\S{t)\\n,,p < \\Sm\n,,p + I (a||<5(s)|U,p + Dr 

and Gronwall's inequality implies 



|<5(t)|U,p<||<5(0)|U,pe'^'-* + ^(e^'-*-l). 



Dr 

It remains to show existence of the above differential equation. This will follow once 
we show that gr+i{t) — gr+i{t) and hr+i{t) — hr+i{t) can be written as a linear 
combination of shifts of 5 with the coefficients depending only on {a{t),b{t)) and 
{a±{t) , b±{t)) . The fact that (a, 5) does not solve TL,. only affects finitely many 
terms and gives rise to the inhomogeneous term Br{t) which is nonzero only for a 
finite number of terms. 

To see that gr+i{t) — gr+i{t) and ~ hr+i{t) can be written as a linear 

combination of shifts of 5 we can use induction on r. It suffices to consider the 
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homogenous case where cj — 0, I < j < r, since aU involved sums are finite. In 
this case |33| Lemma 6.4] shows that gj{n,t), hj(n,t) can be recursively computed 
from gQ{n,t) = 1, ho{n,t) — via 

9j+i{n,t) = ^{h.j{n,t) + hj{n- l,t)) + b{n,t)gj{n,t), 

3 ^ j 

hj+i{n,t) = 2a{n,t)^^gj^i{n,t)gi(n + l,t) - -'^hj^i{n,t)hi{n,t) 

l=Q 1=0 

and similarly for gj{n,t), hj(n,t). Hence the claim follows. 

Finally, observe that since w{n) > 1 this solution is bounded and hence coincides 
with the solution of the Toda equation from Theorem 13.11 □ 



For closely related results we also refer to [35] . 

4. Quasi-periodic finite-gap solutions 

As a preparation for our next section we first need to recall some facts on quasi- 
periodic finite-gap solutions (again see [5], [T3], or [331). 

Let be two quasi-periodic finite-band Jacobi operators^ 

(4.1) H^{t)f{n) = a±(n, t)f{n + 1) + a±(n - 1, t)f{n - 1) + 6±(n, t)f{n) 
in (Z) associated with the Riemann surface of the square root 

23±+l 

(4.2) P±iz) = - n \A^' <Et <---<E^g^^„ 

where G N and ^ is the standard root with branch cut along (— oo, 0). In fact, 

H^{t) are uniquely determined by fixing a Dirichlet divisor J2j=li^^fi^)^'^fi^))^ 
where ^if {t) G -E-l^] aifia^it) G {—1,1}. The time evolution of the Dirichlet 

divisor is determined by the Dubrovin equations (cf. [331 (13.2)] )and linearized by 
the Abel map (cf. [331 Sect. 13.2]). The spectra of H^{t) consist of 5± + 1 bands 

a± 

(4.3) a± a(i/±(t)) ^ [j [Et^,E^^^,]. 

3=0 

We will identify the set C \ a{H^{t)) with the upper sheet of the Riemann surface. 
Associated with H^{t) are the Weyl solutions 

(4.4) ij^iz,n,t)e£\±N) 



normalized such that 'ijj^{z,0,t) — 1. We will use the convention that for A G ct± 
we set -0^ {\,n,t) = limejo "0^ ( A + ie) . Then 

(4.5) i^fi^j t) — exp {af{z, t))tp^{z, rt, t) 
satisfies 

(4.6) 77±(t)^±(z,n,t) = z^±(z,n,t), 

(4.7) ^V'±(z,n,t) = P^^^^^{t)il;^iz,n,t), 



^Everywhere in this paper the sub or super index (resp. "— ") refers to the background on 
the right (resp. left) half-axis. 
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where (03j, (13.47)) 

(4.8) atiz,t) = f (2a±(0,s)G±,(z,0,s)V,±(z,l,s)-ff±,+i(z,0,s))ds. 
Jo 

Note that the integrand in this last expression might have poles if z hes in one of 
the spectral gaps [E^j_i, E^^\. Hence one has to understand a^{z,t) as a limit from 
z e C\R for such values of z. Alternatively one can use the expression in terms of 
Riemann theta functions. Moreover, exp (a^(2;,t)) has simple poles at Ai^(O) and 
simple zeros at fJ.f{t). We refer to the discussion in ^8J for further details. 

5. Inverse scattering transform 

Fix two quasi-periodic finite-gap solutions a^{n,t), b^{n,t) as in the previous 
section. Let a{n, t), b{n, t) be a solution of the Toda hierarchy satisfying 

(5.1) Y.^l + \n\){\a{n,t)-a^{n,t)\ + \b{n,t)-b^{n,t)\) < oo 

n=0 

for one (and hence for any) S I^- In [10] (see also [7], [9], [38]) we have devel- 
oped scattering theory for the Jacobi operator H{t) associated with a{n,t), b{n,t). 
Jost solutions, transmission and reflection coefficients now depend on an additional 
parameter t G M. The essential spectrum of H{t) is (absolutely) continuous and 

(5.2) a(H{t)) = a{H), a,ss{H) = a+iJ a^, (7p(H) = {Afej^^i C R\ae..(H), 
where p e N is finite. We introduce the sets 

(5.3) fj^^^ := (T+ ncr_, cr£^ = clos((T± \ cr(2)), o-:=o-+Ucr_, 

where a is the (absolutely) continuous spectrum of H{t) and a^^^ U cri^'', a^'^^ are 
the parts which are of multiplicity one, two, respectively. 
The Jost solutions ip±{z,n,t) are normalized such that 

(5.4) ^p±{z,n,t) = ^^{z,n,t) {1 + o{l)) as n ^ ±cx). 

Transmission T±{X,t) and reflection R±{X,t) coefficients are defined via the scat- 
tering relations 



(5.5) T^{X,t)ij±{X,n,t) ^ iJj^{X,n,t) + R^{X,t)^P^{X,n,t), A G dzp, 
which implies 



^ _ W{ij±{X,t),^j±iX,t)) ^ _ W{i;^{X,t),i^±iX,t)) 

^ ' ' ^■~mV^T(A,i),V^±(A,i))' ''^^^''>- W{^^{X,t),i;±{X,t)y 

A G a±. Here Wn{f,g) — a[n){f[n)g{n + 1) — f{n + l)g{n)) denotes the usual 
Wronski determinant. 

To define the norming constants we need to remove the poles of i{}^{z,n,t) by 
introducing 

(5.7) ^±(z,n,i) =(5±(z,t)V'±(z,n,t), S±{z,t):= [] {z-^if{t)), 

M,*(t)eM±(t) 

where 

(5.8) M^{t) = {fif{t)\tJ.f{t) GM\a± isapoleofV^(0,l,i)}. 
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The norming constants 7±,fc(i) corresponding to Afe G ap{H) are then given by 

(5.9) 7±,fcW-' =^|^±(Afc,n,i)p. 

nez 

Lemma 5.1. Let {a{t),b{t)) he a solution of the Toda hierarchy such that ()5.1|) 
holds. The functions 

(5.10) V'iC-Z) 1) — exp(a^(z, t))'ip±{z, n, t) 
satisfy 

(5.11) H{t)ip±{z,n,t) ^ z^±{z,n,t), ^ip±{z>n,t) = P2r+2{t)i^±{z,n,t). 

Proof. We proceed as in 31, Theorem 3.2]. The Jost sohitions ip±{z,n,t) are con- 
tinuously differentiable with respect to t by the same arguments as for z (compare 
[71 Theorem 4.2]), and the derivatives are equal to the derivatives of the Baker- 
Akhiezer functions as n — + ±oo. 

For z G C\(T, the solution u±{z, n, t) of (15. lip with initial condition ip±{z, n, 0) G 
^j_(Z) remains square summable near ±oo for all t G R (see [32] or [33l Lemma 
12.16]), that is, u±{z, n, t) — C±{z, t)^(J±{z, n, t). Letting n — !■ ±00 we see C±{z, t) ~ 
1. The general result for all z G C now follows from continuity. □ 

This implies 

Theorem 5.2. Let {a{t),b{t)) he a solution of the Toda hierarchy such that ()5.1|) 
holds. The time evolution for the scattering data is given hy 



T±{\,t) = T±(A,0)exp(a^(A,t) - a?{X,t)), 

(5.12) 



R± (A, t) = R± (A, 0) exp(a± (A, t) - a± (A, t)), 



7±,fc(t) =7±.fc(0)M^^exp(2a±(Afe,t)), I < k < p. 

Proof. The Wronskian of two solutions satisfying (|5.1ip does not depend on n or i 
(see [22, [321 Lemma 12.15]), hence 



T±(A,i) = 



W{i^±iX,t),^±iX,t)) ^ eMafi^t)) W{^P±iX,t),^±iX,t)) 
Wi^^{X,t),^lj±iX,t)) exp(a±(A,0) W{iP^{X,t),i^±{X,t)) 
exp(a^(A,t) - a^iX,t))T±{X,0), 



W{,P^{X),^±{X)) _ exp(a±(A,t))Ty(V;T(A),V^±(A)) 
Wi^^iX),^±iX)) exp(a±(A,0) W{ij^{X),i^±{X)) 
= exp{a^{X,t) - a^{X,t))R±{X,0). 

The time dependence of7±,fc(t) follows from ||C/r(t, 0)?/'±(A/c, ., 0)|| — \\ip±{Xk, .,0)\\. 

□ 



Remark 5.3. Note that we have 
(5.13) 



exp (a± (A, t) - «±(A,t)) = exp ( P± (A) ^f^f '"'J ,, ds ] , A G a± 



n?Si(A-A^f(s)) 
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where Aij^(i) are the Dirichlet eigenvalues. Moreover, in case of the Toda lattice, 
where r — 0, we have G^q (A, n, i) = 1 and H^^ (X,n,t) = A — b^{n,t). 

In summary, since Lemma 13.21 ensures that (|5.ip remains valid for all t once it 
holds for the initial condition, we can compute i?±(A,0) and 7±,fc(0) from (a(n, 0), 
b{n, 0)) and then solve the Gel'fand-Levitan-Marchenko (GLM) equation to obtain 
the sequences (a(n, t),b{n, t)) as in [lOj. More precisely, one needs to solve the GLM 
equation 

(5.14) K±{n,m,t) + Y,K±{n,l,t)F±{l,m,t) ^ ( \v ±m > ±n, 
for /•ir±(n, m, i), where according to Theorem 15 . 21 the kernel F±(m, n, t) is given by 



Theorem 5.4. The time dependence of the kernel of the Gel'fand-Levitan- 
Marchenko equation is given by 

F±{m,n,t) = -Rc [ R± (A, 0)^^^ (A, m, t)^;f (A, n, t) ^'=^ 'u' ^"^^ 



1 ^ ..^nfii(A-M^o))^^ 



(5.15) +— / |r^(A,0)|V±(A,m,i)V'±(A,n,0^^^^i777Yr 

27ri7^(i) ^t(-^) 

+ X! 7±,fc(0)^/'^(Afc, m, i)?/>^(Afc, n, f), 

where i/j^{z, m, t) = (5±(z, 0)';/'^(z, m, i). 

Proof. The kernel F±{m, n, 0) is derived in (TUl Theorem 4.1]. Observe that ajf (A, i) 
are real valued on the set cr^^^ and 



a± 



>^-l^f{t) 



exp(a±(A,t) + a±(A,t)) = ITT^T' ^ ^ a±, 

then our result follows from (|5.10p and Theorem 15.21 □ 

By |10) this equation is uniquely solvable and the solution of the Toda hierarchy 
can be obtained from either Kj^{n, m, t) or K-{n, m, t) by virtue of 



+ ^^ K+{n+l,n+l,t) _ K^{n,n,t) 

ijin,t) = a„ (n,t) ; -. 

«^ ' ^ K+{n,n,t) «^ ' ' K-{n + l,n+l,t)' 



K+{n,n+l,t) _|_ K+{n — l,n,t) 



b{n,t) = b+{n,t)+a+{n,t) ' a+(?i- 1,<) 



K+{n,n,t) « ' ' 'if+(n- l,n- l,t)' 

(5.16) ^ b-in,t) + a-in - M)^^-^;^ - ^ ^ ^) ^_ (J^ ,) • 
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